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Abstract 

We solve the superradiant laser model in two limiting cases. First the sta- 
tionary low-pumping regime is considered where a first-order phase transition 
in the semiclassical solution occurs. This discontinuity is smeared out in the 
quantum regime. Second, we solve the model in the non-stationary regime 
where we find a temporally periodic solution. For a certain parameter range 
well separated pulses may occur. 
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I. INTRODUCTION 



Superfluorescent pulses can be produced by N collectively radiating identical atoms |^-|^ 
as these atoms decay from an initially excited state to the ground state. In contrast to such 
transient behavior would be the stationary output of the superradiant laser discussed in . 
Collective behavior would be manifested in the proportionality of the output intensity to 
A^^ and of the linewidth to A^~^. Moreover, as was shown in such a laser could display 
nearly perfect squeezing of the intensity fluctuations. The theory of a superradiant laser has 
up to now only been considered semiclassically in the stationary regime. As we shall show in 
the present paper, the semiclassical solution needs some quantum mechanical corrections for 
a finite number of atoms and low pumping amplitudes. We also propose to extend previous 
investigations towards a regime not allowing for a time independent stationary solution. 

As in Ref. |^,^ we consider the simplest model of a superradiant laser which accounts 
for A^ three-level atoms placed inside a resonator. Fig. |I]. We assume a resonant coherent 
two-photon pump process 0^2. The two transitions 2 —>■ 1 and 1^0 are coupled to two 
resonant cavity modes a, b, which are assumed to be so strongly damped that they are kept 
in adiabatic slavery by the atoms. 

The atoms are described by the collective population {i = j) and polarisation {i ^ j) 
operators Sij = Yl^=iSij = Z]^(K)(il)^- We shall represent these atomic observables by 
creation and annihilation operators zl,Zi with Sij = zjzj and [2:1,2;]] = 6^; the two modes 
will be described by photon annihilation and creation operators a, b, a^, 6^ The Hamiltonian 
Hq for atoms and field modes reads 

Ho = -ihguiaz^zi - a) z\z2) + i%oi(M^o - 

+ \hn{zlzQ- zlz2). (1) 

Here Vt is the amplitude of the external classical pump field, gi2 and qqi are the atom-field 
coupling constants for the transition 2 1 and 1^0, respectively. As the modes are 
damped we have to add the two irreversible time rates of change for the mode amplitudes 
a, b 




'db' 
dt 



-Kaa{t) + ^^aVait), (2) 



-Kbb{t) + V^bVbit): (3) 



where the two damping constants and T]a{t),r]b{t) the corresponding quantum 

Langevin forces. The latter forces ensure the preservation of the Bose commutators 
[a(t) , a'^ (t)] = [b(t),b^(t)] = 1. The Heisenberg-Langevin equations for the system with 
cavity damping are now 

4 = -gi2azl + goib^zl, (4) 
zl = -nzl - goibzl, (5) 



Z2 



Qzl + guahl, (6) 
a = -gi2ZiZ2 - + y/2K^r]a, (7) 
b = -goizlzi - Kbb + V^bVb- (8) 



2 



These equations have been solved before in under some assumptions that we shortly 
recapitulate here. First the 6-mode is taken to be strongly damped so it can be adiabatically 
eliminated. For later use we define the two damping constants 7a, 7^: 

7a = — , lb = — ■ (9j 

Then an analytic solution can be given in the semiclassical limit ^ 1 in which the 
observables are represented as a sum X = X + 6X of a dominant classical term X ~ and 
a "small" operator valued fiuctuation 6X. The classical term is evaluated in the stationary 
regime by dropping X and the Langevin forces. The resulting solutions can be expressed in 
terms of a dimensionless coupling strength c and an effective pump strength p: 

76 Ny/c-fb 
The solution for the level populations and the field amplitude is 

^ iVc(l-p) ^ 

•^oo = — , Su = Np, 

1 + c 



N{l-p) _ N^, c^p{l-p) 

S22 = ^ , a = , (11) 

1 + c 9i2 Vl + c 

which is stable under the constraints < p < 1 and c > 1 in the case we are interested 
in (additional adiabatic elimination of the a-mode). Notice that for p ^ and p ^ 1 
some of the observables vanish. This is in confiict to the assumption of a large mean value 
and small fiuctuations, so the semiclassical approximation will break down at these points. 
Moreover, this solution cannot be correct in the limit p — >■ since it does not correspond 
to the ground state (5*00) = N. It is physically clear, however, that without pumping all 
atoms will eventually settle in the lowest state, due to damping of the levels 2 and 1. To 
improve on the semiclassical prediction of a discontinuity of (5*22) for p = we have to 
consider the full quantum mechanical solution. This will be done in two ways. First, we 
develop a perturbative small-p expansion. Second, we give a recurrence relation for the 
whole stationary density operator p which can be evaluated numerically. 
In the final section of this paper we shall take a look at the regime where (|1T]) is not a stable 
solution anymore. 



II. FIRST-ORDER PHASE TRANSITION IN THE STATIONARY REGIME 

For the discussion of the quantum mechanical solution we will switch to the master 
equation of the superradiant laser. Additionally we assume the a-mode also to be so strongly 
damped that we can eliminate it adiabatically. The master equation then reads 

p = (L02 + L21 + Lio)p (12) 

with 
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Lo2P = Jl[-So2--520,p], (13) 
L21P = la{[Sl2, PS21] + [5l2P, -521]}, (14) 

LwP = lb{[Soi, pSiq] + [^oiP, -^lo]}- (15) 

For the discussion we expand p in the fully symmetric and normalized states 2') with 
I atoms in level 2, m atoms in level 1, and N — m — I atoms in level 0. The short-hand 
notation |0) = |1°; 2°) is used for the ground state. 

A. Second-Order Perturbation Expansion 

In this subsection wc take a closer look at the neighborhood of the point p = 0. For zero 
pumping the atoms will eventually all settle in the ground state, so the stationary density 
operator is p{p = 0) = |0)(0|. For small non-zero pumping p may be expanded in a series 

p = p(o) + p(i)+p(2) + ... (16) 
where p^") oc p". With A = L21 + -Lio the terms up to second order are 

p(°) = |0)(0|, (17) 



= lim / (itie^(*-*i) Lo2e^*^|0)(0| 



P^'%'){0\ + \0){2'\), (18) 



p(2) = ^limydtie^(*-*^) |, 



t ti 

dt2 


xLo2e^(*i-*^)Lo2e^*^|0)(0| 



+ 

c 



+ (12= X 01 + 10X2^1). (19) 



These entail the mean occupation numbers 



{Soo) ^N- p'N'^ + 0{p'), (20) 
(Sn) ^p'N' + 0(p% (21) 
{S22) ^- + 0{p'). (22) 
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Starting from the ground state the corrections are proportional to and (5*22) is of order 
(as the semiclassical solution) for p = 1/N. Thus the discontinuity of the semi classical 
solution is smeared out over a range < p < 1/N. Further calculations show that the 
fluctuations of the observables vanish for p ^ 0. All of this is reminiscent of a first-order 
phase transition. 



B. Recurrence Relation 

Having studied the perturbation expansion we now prepare to look for a complete quan- 
tum mechanical description. To that end we again expand the density operator p in the 
states |1'";2') introduced above. One may check that coherences with respect to different 
numbers of atoms in level 1 decay to zero. We take this into account by the simple ansatz 

p= a»^,r|l"^;2')(^'";2^. (23) 

0<m + r<JV 
0<m + l<N 

This expansion is inserted into (|12D and we find a recurrence relation for the stationary 

solution PLrn,r 



= pN^/d J{N -m-l + l)lpi-i, 



- ^{N -m-l){l + l)p/+l,m,r 

+ ^J{N -m-r + l)rpi^^^r~i 

- \J {N - m - r){r + l)pi,m,r+i 
+ c |2my'(r + l)(/ + l)p;+ 

l,r?i— l,r+l 

-(m +!)(/ + l)pi,m,r] 



+ 2(m + 1)^{N -m-l){N -m- r)pi,^+i,, 

-m{2N -2m + 2- I -r)pi^rn,r- (24) 

We have not managed to find an analytical solution to this equation, but have solved it 
numerically. Since the number of involved variables is roughly proportional to N^, only up 
to 30 atoms have been considered. The results are shown in Fig. |[ Beside the corrections at 
p = one can make out that the semiclassical solution has to be corrected for p 1 as well, 
as explained above. In the intermediate p range the semiclassical result is quite good, even 
for the moderate values of studied. Since the transition from quantum to semiclassical 
behavior can already be seen in our calculations, it seems not necessary to go to higher N. 



III. NON-STATIONARY REGIME 

In this section we consider the non-stationary solution of our system in the bad-cavity 
limit. As indicated in we expect to find a temporally periodic behavior of the laser field. 
So we set c < 1, where the semiclassical stationary solution is not stable. Starting from the 
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Heisenberg-Langevin equations 
and 2, 



we assume strong, saturated driving between levels 



(25) 
(26) 



^2 = VNsm{nt), 
zq = \/N cos{Qt). 



With the assumption of very fast relaxing cavity modes we can once more eliminate a and 
b adiabatically, 



a(t) 



sm{itt)z{+ W — r]a{t), 

l^a V ^a. 



(27) 



fet(^) = -^5l^sin(^]^)4+ \i^vlit), 



and find a linear equation for zl 



cos(2fit) 



zl + A{t), 



(28) 



(29) 



5 J- 6 



Kb 



(30) 



A(t) = -^2TaSm{nt)r]a{t) + v/2rfe cos(^]t)r/^ (t). 
Now equation (^91) is easily solved 



(31) 



with 



zl{t) = e"(*)zl(0) + e"(*) j e-''^'">A{t')dt' 





^t) = ^^t-^^sin(2ni). 



(32) 



(33) 



Clearly we encounter a periodic behavior with period 1/2Q. We use this as time scale and 
define the dimensionless parameters s, d, and r 



d 



Tb + Tg ^ 1+c 

Vb-Va 1 - c 



2fi 2pv^' 
r = 2fit. 



(34) 

(35) 
(36) 



From this the mean number of photons in the a-mode is 
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(a\t)ait)) = ii^^sin2(r/2) e-'^^-^^'^^M 









Int(r). (37) 



The time integration may be performed after expansion in terms of the Bessel functions 



e 



s sin(r) 



E (-i)«/„(s)e-7 (38) 



and after the death of initial transients e ^ ^ we get 
Int(r) ^ isin2(r/2)e-^^"^(") 

x,E^;^(-.r/„w(i + .-). (39) 

To illustrate the temporal periodicity of this solution we give some examples for various 
parameters p, c in Fig. ^. For the limiting case s ^ 1, which implies small p, we find well 
separated pulses the width and height of which may be estimated in a Gaussian approxima- 
tion. With the parameter e = d/s we find 

1 1 

4^2 



Int(r) = -—-J2np e[2/p-'^(— — - 



with 



cos(r„a^) = -e, sin(r„„^) = -Vl - e^, (41) 

COs{Tmin) = "C, Sm{Tmin) = +Vl - 6^. (42) 

So the pulse width is 



^ /n^^ (43) 

V ^^12^01 

Due to the assumptions of saturated driving and strongly damped cavity modes we have 
^ 3> gi2^9oi^ i^a ^ fl'125 and ^ g^i. This entails the need for a large number of atoms to 
make a small, i.e. to get narrow and well separated pulses. 

The time- dependent spectrum of a light signal which accounts for the width of the filter 
r is given as 0: 

t t 





x(at(ti)a(t2)) (44) 
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Since we have explicit expressions for the a, operators, we plug them in and after some 
integrations sending the width of the filter F — and the observation time t — oo we get 
the following time-averaged expression 



S{u;) 



1603 



EE - 

n I 
In+l+l 

Il+l 



AT Us) (i + f) 



d+'m 



- ) - Un+l 



+ ill 



X 



+ 



d + i{2n + 2l + l + uj/Q) 
1 



d + i(2n + 2/ 



(45) 



We have used the same parameters p, c as in Fig. |^ to illustrate the spectra of some of these 
pulses in Fig. ^. Of course high harmonics of the Rabi frequency VL are enhanced for a ^ 0. 
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FIGURES 



FIG. 1. Scheme of the 3-level superradiant laser. 

FIG. 2. Occupation numbers of the three levels derived through the recurrence relation. Even 
for our moderate number of atoms one can clearly see the transition from the quantum to the 
semiclassical regime. The discontinuity of (-S'22) for p — in the semiclassical solution is smeared 
out over a range IjN. 

FIG. 3. The non-stationary pulsed solution for different parameters. As indicated by the 
calculations the pulse width is given by ^Jp and is independent of c. 

FIG. 4. Spectra of the pulsed solution for different parameters. As one expects the narrow 
pulse {p = 0.1) has higher harmonics than the broad pulse {p = 2). 
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Fig. 1 
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Fig. 2 
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